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Abstract. Hanbury-Brown Twiss (HBT) correlation measurements provide valuable information 
■ about the phase space distribution of matter in ultrarelativistic heavy-ion collisions. The rapidity 

dependence of HBT radii arises from a nontrivial interplay between longitudinal and transverse 
expansion and the time dependence of the freeze-out pattern. For a non-accelerating longitudinal 
expansion the dependence primarily arises from the amount of radiating matter per unit rapidity 
dN/drj, but for a scenario with strong longitudinal acceleration additional complications occur. In 
^sO . this paper I explore schematically what type of dependence can be expected for RHIC conditions 

under different model assumptions for the dynamics of spacetime expansion and freeze-out. 
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INTRODUCTION 

in : 

While the final state of a ultrareativistic heavy-ion collision is to a large degree reflected 
in the measured distribution of hadrons there are different possible evolutions leading to 
Q-i! this final state which leave a more subtle imprint on observables such as transverse mass 
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mj = y m 2 + pj spectra and two particle correlations. 
■ The two extreme assumptions are the boost-invariant model proposed by Bjorken 

Jll] in which the incoming nuclei are nearly transparent and the Landau scenario J2] 
in which the incoming nuclei are completely stopped and the resulting extremely hot 
and dense system is assumed to undergo strongly accelerated expansion. In a Bjorken 
expansion matter moves on free streaming trajectories from a source with negligible 
longitudinal extension (the initial long, nuclear overlap), leading to rapidity r\ = 

being equal to spacetime rapidity r\ s = j In this scenario a plateau in the multiplicity 
per unit rapidity dN/drj is expected around r\ = 0. In contrast, it has been argued that 
dN/drj at RHIC can reasonably well be described with Landau hydrodynamcis O]. In 
this scenario, r] s ^ 77, no plateau is expected and dN/drj changes strongly in time as the 
system undergoes longitudinal acceleration. 

Most 3d hydrodynamic calculations start with initial conditions close to Bjorken- 
like free streaming (albeit without a plateau in rapidity) and consequently do not alter 
dNdrj significantly during the evolution [4]. In contrast, in ||5|] I have suggested that a 
significant longitudinal acceleration seems to be necessary if a simultaneous description 
of transverse mass spectra and HBT correlation radii is required. 

In a longitudinal free-flow scenario (with dN/drj approximately unchanged), the 
effect dominating the rapidity dependence of Hanbury-Brown-Twiss (HBT) radii and 



my-spectra is the shape of dN/drj — the amount of thermalized matter at given T7 
governs the amount of expansion and flow that can be reached before decoupling. 

In contrast, in a scenario in which dN/dr} changes strongly over time, radiation from 
the expanding system will only contribute to the yield at forward rapidities in later 
evolution stages as matter accelerates outward. Thus, e.g. the duration of the presence of 
thermalized matter at forward rapidities will be potentially be smaller than at midrapidity 
quite independent from variations in dN/dr]. However, spectra and HBT radii will only 
reflect this if emission of particles during the evolution time is a significant fraction of 
the total particle emission inlcuding final breakup. It is the aim of this paper to study 
to what degree the rapidity dependence of HBT correlation radii is capable of reflecting 
the evolution history of the system. 



MODEL FRAMEWORK 

I use the model described in for the analysis. Its main assumption is that an 

equilibrated system is formed a short time To after the onset of the collision. This thermal 
fireball subsequently expands isentropically until the mean free path of particles exceeds 
the dimensions of the system and particles decouple at a timescale xy. 
For the entropy density at a given proper time we make the ansatz 

s(x,r ls ,r)=NR(r,x)-H(r, s ,x) (1) 

with T the proper time as measured in a frame co-moving with a given volume element 
and R(r,x),H(r\ s ,x) two functions describing the shape of the distribution and ./V a 
normalization factor. Woods-Saxon distributions 



R(r,x) = 1/ 1+exp 



H{r] s ,x) = l/[ 1+exp 



r\ s -H c (x) 



(2) 



are used to describe the shapes for given x. Thus, the ingredients of the model are the 
skin thickness parameters d ws and 77 ws and the parametrizations of the expansion of 
the spatial scales i? c (r),// c (T) as a function of x. For the radial expansion assuming 
constant acceleration one finds R c (x) = Rq + ^x 2 . H c (x) is obtained by integrating 
forward in x a trajectory originating from the collision center which is characterized by 
a rapidity r\ c (x) = r\Q-\-a-qX with r\ c = atanh v c z where v c z is the longitudinal velocity 
of that trajectory. Since the relation between proper time as measured in the co-moving 
frame and lab time is determined by the rapidity at a given time, the resulting integral 
is in general non-trivial and solved numerically (see [6] for details). Rq is determined 
in overlap calculations using Glauber theory, the initial size of the rapidity interval 
occupied by the fireball matter. The set of model parameters then includes initial rapidity 
T]o transverse velocity = a^Xf and rapidity at decoupling proper time Y\f = riQ+a^Xf. 

Thus, specifying T]o, T]/,Vj_ and Xf sets the scales of the spacetime evolution and J ws and 
T] ws specify the detailed distribution of entropy density. 

For transverse flow a linear relation between radius r and transverse rapidity p = 
atanh v^(t) = r/R c (x) ■ P c (t) is assumed with P c (t) = atanh a±x. The model allows 
for the possibility of accelerated longitudinal expansion which in general implies 77 7^ r} s 



||6|]. This mismatch between spacetime and momentum rapidity can be parametrized as 
a local Arj = rj — rj s which is a function of % and r\ s . 

With the help of a quasiparticle equation of state (EOS), one can find the local temper- 
ature r(7]. s , r, t) from the entropy density s{r\ s , r, r) and net baryon density Pb(Vs, r, t) . 
Particle emission is calculated throughout the lifetime of the fireball by selecting a 
freeze-out temperature Tf, finding the hypersurface characterized by T(r} s ,r,l) = Tf 
and evaluating the Cooper-Frye formula 



d 3 N 



d 3 p (2n) 



do^p^ exp 



d 4 xS(x,p) 



(3) 



with p^ the momentum of the emitted particle and g its degeneracy factor. 

HBT correlation radii are determined as averages of coordinates with the emission 
function 0] 



R 2 side (K) = (f)(K), R 2 out (K) = {(x-PJ) 2 )(K) and Rj ong (K) = (z 2 )(K) (4) 
where 

By choosing rjo different from the measured width of dN/drj it is possible to intro- 
duce longitudinal acceleration into the model. Note that the model description cannot 
be assumed to be valid in the target/projectile fragmentation — as one probes forward 
rapidities, thermal physics is less relevant, evident from the fact that a large fraction of 
matter is below Tp ab initio. Thus, the model prediction beyond r\ = 2 should be taken 
with some caution. 



Rout AND THE FREEZE-OUT PATTERN 

We expect that the spacetime evolution of the emission source is reflected in the rj 
dependence of HBT radii only if the system is not dominated by the final breakup. Since 
freeze-out is determined dynamically, there is no apparent parameter to adjust amount of 
hadrons released by the final breakup vs. the amount of particles continuously emitted 
before. From Fig. [T] it is however evident that selecting the skin thickness d ws in the 
radial distribution of entropy density has a visible impact on the radial position Rcf(^) 
of the Cooper-Frye surface. 

Rcf shrinks for a large skin thickness (almost a Gaussian) whereas it grows for a 
small skin thickness (a box-like distribution). Since particles must move faster than the 
Cooper-Frye surface in order to be released, an inward-moving surface implies much 
more early emission than an outward moving one. This is confirmed in the right panel 
of Fig. [T] where we present the number of particles being released as a function of x. 
The width of the final peak in this figure is connected with the emission time, i.e. the 
difference in R (mt and R s Me- Note that the data seem to prefer the solution with expanding 
Rcf, i-e. a sharp surface, at least in the last stages of the evolution. Shockwaves in the 
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FIGURE 1. Left panel: Evolution of the Cooper-Frye surface radius R C f at tj — for 7> = 1 10 MeV 
as a function of T for different values of radial skin thickness parameter. Right panel: Particle emission (in 
arbitrary units) into midrapidity as a function of T for pions as a function of Rcp{l) and for photons. 
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FIGURE 2. R s ide as a function of pair momentum kj for different rapidities rj in an approximate 
boost invariant expansion (left), a longitudinally accelerating scenario with continuous particle emission 
(middle) and a longitudinally accelerating scenario with sudden breakup (right). 



medium excited by hard parton energy loss as observed in two particle correlations 
could be a potential mechanism for creating such a sharply defined surface. 



RESULTS 



I investigate rapidity dependence in three different scenarios designed to show the 
essential effects to be expected: 1) a longitudinal free-flow expansion where dN/dr} 
remains approximately constant during the evolution 2) a longitudinally accelerated 
scenario where the width of dN/dr} changes from 1.7 to 3.8 during the evolution and 
continuous particle emission occurs and 3) a longitudinally accelerated expansion with 
the same parameters as above but a sudden breakup in the end (note that case 1) is also 
characterized by a sudden breakup). In all three cases single particle mj spectra describe 
the data at midrapidity and the measured dN/dr] is reproduced, though only 3) describes 
the measured HBT correlation radii at midrapidity correctly. 

Figure |2l shows R s ide in the three different scenarios as a function of pair transverse 
momentum for different rapidities r) . 




1- 1- - 1- 



$.2 0.4 0.6 0.8 1 8.2 0.4 0.6 0.8 1 8.2 0.4 0.6 0.8 

kj, [GeV] k,. [GeV] kj, [GeV] 

FIGURE 3. R out as a function of pair momentum kj for different rapidities r\ in an approximate 
boost invariant expansion (left), a longitudinally accelerating scenario with continuous particle emission 
(middle) and a longitudinally accelerating scenario with sudden breakup (right). 



The differences between the three scenarios is generally small, however the geomet- 
rical radius in the free-flow expansion is systematically smaller than in the accelerated 
cases. This is due to the fact that breakup occurs roughly at a fixed volume (determined 
by the total entropy), hence the stronger longitudinal expansion in the Bjorken case 
(which expands with T]o at all times) implies less transverse radius expansion and flow. 

Fig. E] compares R out for the three different scenarios. The most striking effect is the 
general increase in the free flow and the continuous emission scenario as compared to 
the sudden breakup one. The reason can easily be identified: With an expanding Cooper- 
Frye surface and a large fraction of energy going into transverse expansion, the sudden 
breakup scenario has the strongest positive x — t correlation and the least release of 
particles prior to final breakup. While sudden breakup also occurs in the approximate 
Bjorken solution, the stronger longitudinal expansion leaves less energy for transverse 
expansion, significantly lessening this correlation. Finally, since the Cooper-Frye surface 
is moving inward in the continuous emission scenario, there is no such correlation and 
R out significantly exceeds R side . 

Finally, in Fig. HJl compare the three different scenarios with respect to Ri ong . Due 
to r\ « rj s in this scenario but r\ > ri s in the presence of longitudinal acceleration, the 
scaled falloff with kj is the same in all three cases (since the final dN/drj is fixed), 
the correlation radius normalization however is given by r} s and this reduces the overall 
normalization. The radii in the sudden breakup scenario are systematically larger than in 
the continuous emission one, reflecting the fact that freeze-out happens on average later 
and hence at larger spatial extension. The sudden drop at r\ = 3 is probably an artefact 
of the model getting sensitive to the target/projectile fragmentation region which is not 
well described by thermal physics. 

CONCLUSIONS 

The most remarkable finding is that to first order there is hardly any influence of the 
expansion and freeze-out scenario on the rapidity dependence of HBT correlations. The 
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FIGURE 4. /?/ on g as a function of pair momentum kj for different rapidities rj in an approximate 
boost invariant expansion (left), a longitudinally accelerating scenario with continuous particle emission 
(middle) and a longitudinally accelerating scenario with sudden breakup (right). 

reason for this is that the measured single particle spectra and dN/drj already provide 
rather stringent constraints on the expansion pattern. The most pronounced effect of a 
longitudinally accelerated expansion is the deviation from the scaling r\ m r} s , ultimately 
leading to a smaller Ri on g m me presence of acceleration while keeping the correct falloff 
in kj- For the most part changes in the sideward and outward correlation radii are minor. 
Since the data practially rule out strong continuous emission throughout the lifetime 
by Rout/Rside ~ 1> the correlation radii are dominated only by the state of the system 
at breakup, not by the evolution path. It seems that for reasonable evolution scenarios 
which are in agreement with other hadronic data the r\ dependence of HBT radii is 
primarily determined by the amount of matter at given r\. Thus, photons might be a 
more promising tool in order to study differences in the evolution history [10]. 
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